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^ ■ ABSTRACT 

bJQ, Context, Double-periodic (beat) Cepheids are important astrophysical objects which allow testing both the stellar evolution and 

stellar pulsation theories, as well as the physical properties of matter in stellar conditions. However, the phenomenon of double- 
periodic pulsation is still poorly understood. Recently we rediscussed the problem of modelling the double-periodic pulsation with 
non-linear hydrocodes. We showed that the published non-resonant double-mode models are incorrect, because they exclude the 
C*") , negative buoyancy effects. 

Aims. We continue our efforts to verify whether the KuhfuB one-equation convection model with negative buoyancy included can 
reproduce the double-periodic Cepheid pulsation. 

Methods. Using the direct time integration hydrocode, which implements the KuhfuB convection model, we search for stable double- 
periodic Cepheid models. We search for models pulsating in both fundamental and first overtone modes (F+IO), as well as in the 
two lowest order overtones (10+20). In the latter case, we focus on reproducing double-overtone Cepheids of the Large Magellanic 
Cloud (LMC). 

Results. We have found full amplitude non-linear beat Cepheid models of both types, F+IO and 10+20. In the case of F+IO 
models, the beat pulsation is most likely caused by the three-mode resonance, 2ll>i = loq + ati, while in the double-overtone models 
the underlying mechanism (resonant or non-resonant) cannot be identified beyond doubt. Double-periodic models found in our survey 
exist, however, only in narrow period ranges and cannot explain the majority of the observed double-periodic objects. 
' Conclusions. With only little doubt left, we conclude that current one dimensional one-equation convection models are incapable 

of reproducing the majority of the observed beat Cepheids. Among the shortcomings of current pulsation hydrocodes, the simple 
treatment of convection seems to be the most severe one. Growing evidence for the presence of non-radial modes in Cepheids suggests 
that the interaction between radial and non-radial modes should also be investigated. 

■ Key words, hydrodynamics - convection - methods: numerical - stars: oscillations - stars: variables: Cepheids 
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^) 1 . Introduction during the blue-ward and then the red-ward evolution. The cor- 

; . responding second and third crossings last roughly two orders 

OO Cepheids are periodic variable stars, pulsating radially with pe- of magn it u de longer than the first crossing. Hence, most of the 

O riods ranging from a single day to over one hundred days for obser ved Cepheids are expected to be helium-burning objects. 
t^> the most luminous variables. They obey an empirical period- 

. luminosity relation, which ranks them among the most impor- Most of the Cepheids are single-periodic variables pulsating 

tant distance indicators in astrophysics. They populate a rather either in the fundamental mode or in the first overtone mode, 

narrow strip in the HR diagram, called the classical instability Only a few examples of second overtone Cepheids are known, 

strip, in which pulsation is driven by the opacity mechanism act- In many Cepheids simultaneous pulsations in two modes are ob- 

H \ ing in hydrogen-helium partial ionisation region. They are giant, served. Double-periodic Cepheids (or beat Cepheids) are very 

relatively cool stars with typical effective temperatures ranging interesting and important astrophysical objects. Most of these 

from log r e ff = 3.7 to log r eff = 3.8 (see e.g. Sandage, Tammann variables pulsate simultaneously in two consecutive radial pul- 

& Reindl 2004). Estimated values of the width of the instability sation modes, either in the fundamental and in the first overtone 

strip do not exceed 1000 Kelvins Q Depending on the mass, the (F+IO) or in the two lowest order overtone modes (lO-t-20, 

star can cross the instability region up to three times. The first double-overtone Cepheids in the following). The two pulsa- 

crossing occurs in the post-main sequence evolution, before he- tion periods can be used to constrain the stellar parameters; 

lium ignites, while the star evolves very quickly towards the red in particular, the mass can be derived using the Petersen di- 

giant branch. After the helium ignition, more massive stars enter agram (Petersen 119731 1. The disagreement between these pul- 

the horizontal blue loop, crossing the instability strip twice, first sation masses of F+IO beat Cepheids and the masses derived 

from evolutionary computations was one of the factors that mo- 
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The transformation from the observed colours to effective temper- tlvated the revision of opacity tables (SimonQjjSl. Multi-mode 

atures is particularly difficult for pulsating stars. In addition, inaccurate Cepheids still provide useful and stringent tests for stellar evo- 

or unknown reddenings for individual objects significantly add to the lution and pulsation theories (see e.g. Moskalik & Dziembowski 

uncertainty of the log T cS calibration. 120051 Dziembowski & Smolec 2009 ). 
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Despite the great importance of the beat Cepheids and the ap- 
parent simplicity of their oscillation, the phenomenon of double- 
periodic pulsation is not well understood. Many efforts were 
made over the past forty years to study the interaction of pul- 
sation modes and the origin of the double-periodic pulsation. 
These include theoretical studies as well as numerical modelling. 

The necessary condition for the double-periodic pulsation 
to occur is simultaneous linear instability of the respective two 
modes. This condition is not sufficient because mode selection 
is a non-linear phenomenon. The interaction between pulsation 
modes may be either of a non-resonant or of a resonant char- 
acter. In most of the observed Cepheids both mechanisms lead 
to stable single-periodic pulsation. In the non-resonant case, it 
is also a single-mode pulsation. The distinction between multi- 
periodic and multi-mode pulsation is important. In the resonant 
case the single-periodic pulsation is actually multi-mode, be- 
cause the resonantly coupled mode is also excited to high am- 
plitude. Owing to the non-linear frequency synchronisation, it 
does not appear with a separate frequency, but its presence is 
manifested in the distortion of light and radial velocity curves. 
Bump Cepheids, which are single-periodic, double-mode pul- 
sators, provide an excellent example. In this paper we deal with 
double-periodic pulsators and will use the term double-mode be- 
low only for stars and models for which the resonant mechanism 
is excluded. 

Double-periodic Cepheids are rare compared to single- 
periodic Cepheids, and the identification of the conditions that 
lead to stable double-periodic pulsation was and still is of key 
importance. Simon (119791 1 suggested that the double-periodic 
pulsation may be a resonant phenomenon, connected with the 
three-mode resonance, a>i + a>o = ojt,. Using the amplitude 
equation formalism Dziembowski & Kovacs ( 1984) showed that 
the resonance proposed by Simon actually stabilizes the single- 
periodic pulsation. They also pointed out that the 2:1 resonance 
between the linearly excited mode and the linearly damped, 
parasite mode can lead to double-periodic pulsation. This was 
later confirmed with radiative RR Lyrae and Cepheid models 
(Kovacs & Buc hler [19881 Buchler, Moskalik & Kovacs [19901 
Smolec 2009a). The models agreed not once with the obser- 
vations, though. Indeed, a comparison of the observed periods 
and the period ratios of beat Cepheids with the results of linear 
modelling indicate that in most of the double-periodic Cepheids 
resonances of a low order are excluded, and therefore the non- 
resonant mechanism should be operational. For many years ra- 
diative models failed to reproduce the non-resonant double- 
mode pulsation in both RR Lyrae stars and in Cepheids. It was 
the inclusion of turbulent convection into the hydrocodes that 
led to success. Kollath et al. (1998) published the double-mode 
Cepheid models and Feuchtinger ( 1 19981 ) published a double- 
mode RR Lyrae model. Later, only the Florida-Budapest group 
computed the surveys of double-mode pulsation (Kollath et al. 
120021 Szabo, Kollath & Buchler [20041 Buchler 120091 ). The re- 
liability of these models was recently questioned however by 
the analysis of Smolec & Moskalik d2008al [2008b) (see also 
Sect. 12.21 ). We have shown that the computed double-mode 
pulsation arises from the neglect of negative buoyancy in the 
Florida-Budapest hydrocode, which is physically not justified. 
Yet we could not find any non-resonant F+10 Cepheid mod- 
els even with our pulsation hydrocode (Smolec & Moskalik 
2008a, Sect. 12.11 ). which properly included the negative buoy- 
ancy (Smolec & Moskalik l20U8b1) . 

In the present paper we discuss some resonant F+IO 
Cepheid models that we computed (Sect. [3). Although these 
models are restricted to a very narrow parameter range, they 



agree sufficiently well with the observations. Finally we present 
a survey of models which we conducted in search for stable 
double-overtone (10+20) Cepheid pulsation (Sect. [4]). This sur- 
vey is complementary to our earlier work, in which we searched 
for non-resonant F+IO Cepheid models (Smolec & Moskalik 
2008b), and allows us to draw more definite conclusions about 
the ability of current convection models to reproduce the beat 
Cepheid pulsation (Sect. [5]). Preliminary results of this work 
were presented in Smolec & Moskalik (2009). 

2. Turbulent convection model and numerical 
methods 

2.1. Turbulent convection model and pulsation hydrocodes 

In all our computations we use the pulsation codes described by 
Smolec & Moskalik (2008a). These are a static model builder, 
linear non-adiabatic code and a direct time-integration non- 
linear hydrocode. The codes use a simple Lagrangian mesh. 
For the convective energy transfer we use the time-dependent 
KuhfuB C 1986) convection model reformulated for the use in 
stellar pulsation codes (see also Wuchterl & Feuchtinger 1998 ). 
Radiation is described in the diffusion approximation. Below we 
provide a short summary of the model. For extensive description 
and details of numerical implementation we refer the reader to 
Smolec & Moskalik d2008al) . 

The momentum, internal energy, and turbulent energy equa- 
tions are 



du Id, x GM r 

17 = —7r\P + P i ) + u ^ — r 

dt p or v ; r L 
dE dV _ ld[r 2 (F I + F c )] 



— + p — 

dt 1 dt 
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r 2 dr 



C. 



de, 



dV 



d7 + A d7 



1 d(r 2 F t ) 
p r 2 dr 



+ E q + C. 



(1) 



(2) 



(3) 



Above, u is the fluid velocity, which is a time derivative of the 
radius u = dr/dt. M r is the mass enclosed in the radius r, V is 
the specific volume which is the inverse of the specific density, 
V = I /p. p and E are the pressure and energy of the gas. F r , F c 
and F t are the radiative, convective, and turbulent fluxes, respec- 
tively. The radiative flux is computed assuming diffusion approx- 
imation. Radiation pressure and radiation energy are included in 
p and E. The turbulent energy, e t , is computed according to the 
one-equation model of KuhfuB (Eq. 0. Turbulent energy equa- 
tion and internal energy equation are coupled through the term 
C of the form 



C = S — D — D t , 
where, 

c TpQ L/2 

S = aa s ——e t ' Y, 
H p 



D = a d 



j_t 

aH p ' 



4o- Tr 2 r 3 v 2 

Dr = — 2 m et - 



(4) 



(5) 



(6) 



(7) 



Above, T is the temperature, cr is the Stefan-Boltzmann con- 
stant, Q = (dV/dT) p , H p is the pressure scale height, c p is the 
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Table 1. Three sets of convective parameters considered in this 
paper. Parameters a s , a c , a^, Q'p and y r are given in the units of 
standard values (a s = a c = 1/2 V2/3, = 8/3 y/2/3, a p = 2/3 
and y r = 2 V3; see text below and Smolec & Moskalik 2008~alfor 
details). 



Set 


a 


ce m 




a c 


ttd 


a p 


a t 


7r 


Rl 


1.5 


0.30 


1.0 


1.0 


1.0 


1.0 


0.0 


0.75 


R2 


1.5 


0.30 


1.0 


1.0 


1.0 


0.5 


0.0 


1.00 


R3 


1.5 


0.50 


1.0 


1.0 


1.0 


0.0 


0.0 


1.00 



2.2. Importance of negative buoyancy 

In this section we look in detail at the important difference be- 
tween the convection model adopted in our code and the con- 
vection model adopted in the Florida-Budapest hydrocode (e.g. 
Kollath et al. 2002 ). The latter code was used to compute all the 
convective double-mode Cepheid models published so far. As 
analysed by Smolec & Moskalik (2008b), these double-mode 
models are physically incorrect, because they result from the ex- 
clusion of the buoyant forces in convectively stable regions. In 
the Florida-Budapest hydrocode it is set 



specific heat at a constant pressure and k is the opacity. The 
source (or driving) function, S , describes the rate of turbulent en- 
ergy generation/damping through the buoyant forces. The source 
function is proportional to the superadiabatic gradient (dimen- 
sionless entropy gradient), Y, 



Y = V - V a 



H^ch 
c p dr' 



(8) 



and therefore drives the turbulent energies in convectively un- 
stable regions (S > 0), and brakes the turbulent motions in con- 
vectively stable regions (S < 0). Term D models the decay of 
turbulent energy through the turbulent cascade. D r describes the 
rate at which the turbulent energy is transformed into the internal 
energy through the radiative cooling of the eddies (see Wuchterl 
& Feuchtinger 1998). For the turbulent fluxes we have 



F c = aa c pTc p e^ 2 Y, 

Fi = -aa t pH p e t ' — , 
dr 



p t = a v pe t , 



U 



J_d_ 

pr^ dr 



-aa m pH p e\ l2 r ' 
3 \ or r 



du 



E q = -aa m H p e\ l2 \ - - 
3 \or r 



du u 



(9) 
(10) 
(11) 
(12) 

(13) 



Above, pi is the turbulent pressure and U q and E q are the viscous 
momentum and energy transfer rates. We note that the turbulent 
viscosity always contributes to the driving of the turbulent en- 
ergy, at the cost of pulsation. 

The model equations contain 8 orders of unity scaling pa- 
rameters, mixing-length, a and parameters multiplying the tur- 
bulent fluxes and driving/damping terms, a p , a m , a c , a u a s , 
and y r . The theory provides no guidance for their values, but 
some standard values are in use. These values result from a com- 
parison of the static, time-independent version of the model with 
the standard mixing-length theory (see Wuchterl & Feuchtinger 
1998 and Smolec & Moskalik l2008al l. In practice, the values of 
the parameters should be determined in a way that the models 
satisfy as many observational constraints as possible. In prin- 
ciple, direct numerical simulations of convective zones can be 
used to calibrate the convective parameters, but currently no such 
simulations are available for large amplitude giant pulsators. The 
values of the parameters used in the present paper are given in 
TableQ] We note that set R3 is equivalent to set B of Baranowski 
et al. (120091 1 and was adopted by these authors in successful 
modelling of the overall properties of the radial velocity curves 
of first overtone Galactic Cepheids. 



S cc y + , 



for the turbulent source function, and 



F,ocY + 



(14) 



(15) 



for the convective heat flux, to be compared with S cc Y and 
F c cc Y in our code (Eqs.[5]and[9]). To exclude the buoyant forces 
in convectively stable regions, a consequence of Eq. ( TBI ) in the 
Florida-Budapest model, is incorrect and cannot be justified. In 
convectively unstable layers (Y > 0) the buoyant force is respon- 
sible for the acceleration of turbulent eddies. In convectively sta- 
ble regions (Y < 0), the buoyant force does not disappear. It de- 
celerates the turbulent eddies, restoring the convective stability. 
Thus, we stress that the treatment of the source function in the 
convectively stable regions, S cc Y in our code and S cc Y+ in the 
Florida-Budapest code, is not arbitrary. Even more importantly, 
the use of two different equations for the source function (Eq. [5] 
or[T4l leads to significant differences in the computed models. 
Particularly, the modal selection is very different. As Smolec & 
Moskalik (2008b ) showed, when negative buoyancy is neglected 
(Eq. fT4l the non-resonant F+IO double-mode models can be 
easily computed. But the mechanism responsible for the double- 
mode pulsation is artificial. Because the negative buoyancy is ex- 
cluded, the turbulent eddies cannot be braked effectively below 
the envelope convection zone. Consequently, a large region ex- 
tending over several pressure scale heights below the convection 
zone with relatively high turbulent energies is present in these 
models. In this region eddy-viscous damping strongly limits the 
amplitude of the fundamental mode; much stronger than it limits 
the amplitude of the first overtone. With its reduced amplitude, 
the fundamental mode is no longer able to saturate the pulsation 
instability alone, which allows the growth of the first overtone 
and leads to stable double-mode pulsation. For more details the 
reader is referred to Smolec & Moskalik (2008b) and Smolec 
(I2009bt . 

Our extensive survey (Smolec & Moskalik 2008b ) in search 
for stable F+ lO double-mode pulsation in the physically correct 
models, i.e. with negative buoyancy included, yielded no result. 
We could not find a non-resonant double-mode model. However, 
we have found some interesting F+IO resonant double-periodic 
models, which we describe in Sect. [3] We also performed 
the model survey intended to reproduce the observed double- 
overtone Cepheids in the LMC. Our results, also including 
double-overtone models, are described in Sect. [4] 



2.3. Mode selection analysis 

To search for stable double-periodic pulsation we adopt the 
methods developed by the Florida-Budapest group (e.g. Kollath 
et al. 120021 1. They are based on the analysis of results of direct 
hydrodynamical model integration with an amplitude equation 
formalism (e.g. Buchler & Goupil 1984). Here we provide a brief 
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summary of the method. For details of our implementation the 
reader is referred to Smolec & Moskalik (2008b). 

The static model is initialized (or "kicked") with a scaled 
mixture of linear velocity eigenvectors of the studied modes. 
Then, a time-evolution of mode amplitudes follows using the an- 
alytical signal method (see e.g. Kollath et al. 2002). The model 
integration is stopped when full amplitude single-periodic pul- 
sation (limit cycle) is reached, or when the time evolution be- 
comes prohibitively slow, which may indicate the approach to 
the attractor. The resulting trajectory is plotted in the amplitude- 
amplitude diagram, similar to those shown in Fig.Q] In this fig- 
ure we plot on the vertical axis the fractional radius amplitude, 
6R/R, of the first overtone, and on the horizontal axis we plot the 
fractional radius amplitude of the fundamental mode. Each tra- 
jectory is represented by a single, continuous line starting in the 
bottom left part of the diagram, with the direction of the evolu- 
tion marked by an arrow. For each model the non-linear integra- 
tion is repeated with several different initial conditions (range 
of amplitude ratios), chosen to adequately explore the whole 
amplitude-amplitude phase-space. The dimension of the phase- 
space can be large, depending on the number of linearly excited 
modes and possible resonances. Linearly damped modes that can 
be present in the initial phases of the model integration (due to 
imperfect initialization) decay exponentially and are of no in- 
terest, unless they are in resonance with other linearly excited 
modes. In the latter case, their amplitude can grow significantly. 

Provided that the amplitude-amplitude phase-space is suf- 
ficiently covered (several trajectories, sufficiently long integra- 
tions), the mode selection can be deduced beyond doubt from the 
analysis of the appropriate 2D amplitude-amplitude diagrams 
similar to those shown in Figs. Q] and [6] (which are the sections 
through the whole phase-space). For example, the trajectories 
presented in first panel of Fig. UKr^ = 5850 K) evolve either 
towards fundamental mode single-periodic pulsation (a funda- 
mental mode limit cycle; two rightmost trajectories) or towards 
first overtone single-periodic pulsation (six leftmost trajectories, 
one of them overlaps with the vertical axis). In four consecutive 
panels of Fig. [1] the double-periodic attractor is clearly visible. It 
is also evident that for these models the fundamental mode limit 
cycle is stable, while the first overtone pulsation is unstable. For 
the last model presented in Fig.Q](r e ff = 5800 K), all trajecto- 
ries evolve towards fundamental mode, which is the only stable 
pulsation state. More robust conclusions can be drawn however 
if the hydrodynamic trajectories are analysed with appropriate 
amplitude equations. 

Amplitude equations can be used to compute the time- 
evolution of the mode amplitudes and phases. In particular, the 
mode selection for a given model can be established through 
computing the time-independent solutions of the amplitude 
equations (fixed points) and their stability. Stable fixed points 
are the attractors of the system and the trajectories evolve toward 
them. Unstable fixed points repel the trajectories. Fixed points 
correspond to steady non-linear pulsations. Stable, single-mode 
fixed points correspond to stable limit cycle pulsations. 

In order to compute the fixed points and their stability, 
the linear growth rates as well as non-linear self- and cross- 
saturation coefficients of the modes have to be known. In a reso- 
nant case it is also necessary to know the coupling coefficients of 
the interacting modes. These coefficients can be computed from 
the linear eigenvectors and model structure, but this is too diffi- 
cult, unless strongly simplifying assumptions are made. In prac- 
tice, saturation and coupling coefficients can be derived through 
fitting the appropriate amplitude equations to non-linear trajec- 
tories computed with the hydrocode. The numerical procedure 



is straightforward in a non-resonant case. The amplitude equa- 
tions, which are complex in general, decouple into a real part 
that yields the equations for the amplitudes of the modes, and 
an imaginary part that yields the equations for their phases. The 
equations for the amplitudes are used to derive the saturation co- 
efficients through a simple linear fit, while the equations for the 
phases are not relevant in the non-resonant case. 

In a resonant case, the situation is much more complicated. 
Full complex amplitude equations have to be considered. The 
number of unknown coefficients grows significantly, and a sim- 
ple linear fit is no longer possible. Therefore, the mode selec- 
tion analysis, as described above, is hard to conduct (see Smolec 
2009b). Below the modal selection for models in proximity of 
the resonances is derived based on the analysis of hydrodynamic 
trajectories only. As mentioned above and presented in the forth- 
coming sections, the computation of several trajectories and suf- 
ficiently long integrations allows us to deduct the modal selec- 
tion beyond doubt also in the resonant case. 



3. Resonant beat F+10 Galactic Cepheid models 

There are only a few double-periodic Cepheids known in our 
Galaxy. The high extinction in the direction of the Galactic 
disc, where we expect to find these young objects, signifi- 
cantly limits our detection ability. Including very recent discov- 
eries in all-sky surveys data, such as AS AS (Poimanski I2002I ). 
there are 23 known F+IO double-periodic Cepheids (Antipin 
IT9971 [19981 l2006l Ber dnikov & T urner [19981 Khruslov 1200911 
Pardo & Poretti [19971 Wils et al. l20T0l Wils & Otero [20041 . 
Only 15 10+20 Cepheids were detected so far (Beltrame & 
Poretti 120021 Hajduet al. [20091 Khrus lov I2009bl I2009cl [20101 
Pardo & Poretti [19971 Szczv del 120091 ). of which only two ob- 
jects were known before 2009. In many cases the reddenings 
are unknown, making the derivation of physical parameters of 
these objects very uncertain. Here we focus on F+IO Cepheids. 
Observed period ratios, P\/Pq, are confined between 0.697 and 
0.713 for most of the objects. The typical periods of the fun- 
damental mode are between two and six days (see Moskalik & 
Kolaczkowski 2009). The period ratio decreases with increasing 
period of the fundamental mode. 

The double-periodic F+IO models presented in this section 
were found accidentally during the model survey intended to 
reproduce the Hertzsprung bump progression (Smolec 2009b). 
A detailed linear analysis revealed that the computed double- 
periodic behaviour can be connected to the high-order paramet- 
ric resonance, 2a>i = a>o + ol>2- The same resonance was found 
to be responsible for the double-periodic behaviour computed 
in radiative f3 Cephei models by Smolec & Moskalik (120071 ) . 
/3 Cephei stars are multi-periodic, mostly non-radial pulsators, 
and two radial modes excited alone were so far not observed 
in any star of this class. Consequently, the double-periodic 
/3 Cephei models of Smolec & Moskalik (2007 ) are of theoretical 
interest only. Yet they demonstrate that the three-mode, high- 
order resonance can be conducive to producing stable multi- 
mode pulsation. The convective 6 Cephei models described be- 
low confirm this finding and agree moreover fairly well with the 
observations. 

Static models were constructed using 150 mass zones, ex- 
tending down to 2.5 x 10 6 K, with a fixed temperature (T a = 
1 1000 K) in the anchor zone located 50 zones below the surface. 
As discussed in Smolec & Moskalik d2008al ). T a = 15000 K was 
used for model sequences in which effects of turbulent pressure 
were turned on (a p ^0). The convective parameters of the model 
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Fig. 1. Fractional radius amplitude for the first overtone (A{) and for the fundamental mode (Aq) for six consecutive Cepheid 
models of set R2. The models were initialized with a range of amplitude ratios (bottom left parts of the diagrams) and for each 
initialization the time evolution (marked with arrow) was followed. The model properties are given in each panel. The proximity 
to the resonance centre is characterised by A = 2wi/(«o + oji)- Solid and open squares mark the location of stable and unstable 
fixed points computed through fitting the hydrodynamical trajectories with non-resonant amplitude equations. In panels three to five 
some trajectories evolve to a double-periodic attractor with finite amplitude in both the fundamental and first overtone modes. The 
remaining (rightmost) trajectories calculated for these three models evolve towards a single-periodic finite amplitude fundamental 
mode pulsation. 



sequences that are discussed in this section, Rl, R2 and R3, are 
collected in Table Q] 

The computed models are characterised by Galactic chem- 
ical composition, X = 0.7, Z = 0.02. OPAL opacities were used 
(Iglesias & Rogers 1996) and were supplemented at the low tem- 
peratures with the Alexander & Ferguson (1994) opacity data. 
Opacities were generated for the solar mixture of Grevesse & 
Noels ( 119931 1. Models were computed along a sequence of con- 
stant mass and constant luminosity, and varying effective tem- 
perature. Masses of the models are 4.5 M Q for models adopt- 
ing convective parameters of set Rl, 4.75 M Q for models adopt- 
ing convective parameters of set R2, and 5.0 M Q for models 
adopting convective parameters of set R3. Luminosities of the 
models were derived using the mass-luminosity relation result- 
ing from the Schaller et al. (1992) evolutionary computations, 
log(L/L ) = 3.56 log(M/M G ) + 0.79. 

The multi-mode solution was first found for one model of 
set Rl. Detailed linear analysis revealed that the 2cj\ = wo + a>2 
resonance can be involved in the pulsations, because the model 
was located very close to the resonance centre. The proximity 
parameter, which we define for the resonance under discussion 
as 



was equal to A = 1.0006. This finding motivated the search 
for double-periodic F+IO Cepheid pulsation connected with 
the discussed three-mode resonance. Multi-mode models were 
found also for sets R2 and R3 of Table Q] The widest domain of 
double-periodic solutions exists for parameter set R2. In Fig. Q] 



Table 2. Properties of the computed double -periodic F+IO con- 
vective Cepheid models (mass, effective temperature, proximity 
parameter, fundamental mode period and P\/Pq period ratio). 
Convective parameters of the models are given in Table Q] 



Set 


M[M ] 


7df[K] 


A P [d] 


Pl/Po 


Rl 


4.50 


5755 


1.00064 4.21 


0.6964 


R2 


4.75 




see Fig. [T] 




R3 


5.00 


5965 


1.00009 4.77 


0.6929 



we present the results of hydrodynamical model integrations for 
six consecutive models of this set. In Fig. [2] the loci of the 
2a>i = o)()+(l>2 resonance in the HR diagram are plotted. Double- 
periodic models are located within the circle, at log(L/L Q ) « 3.2. 
In Fig. |2] we also plot the lines representing the loci of other res- 
onances, which are important for Cepheid pulsation, 2wo = C02 
and 2u>\ — L04. 

The mode selection scenario is evident from the hydrody- 
namical trajectories alone and is the same for all three parameter 
sets under discussion (see Fig. [T). The double-periodic attrac- 
tor always coexists with the stable fundamental mode attractor. 
The hysteresis domain (in which the mode selection depends 
on initial conditions) is always very narrow. Close to the reso- 
nance centre the consecutive models were computed in 5 K steps 
in the effective temperature. For sets Rl and R3 the double- 
periodic solution was found in only one model (see Table [2] 
for its properties), indicating that the interesting domain is nar- 
rower than 10 K. Only for models adopting convective parame- 
ters of set R2 the double-periodic domain is wider, and extends 
for more than 40 K. For all three model sequences, this domain 
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Fig. 2. Location of the 2w\ = cjq + u>2 resonance centre (dot- 
ted line) and the computed resonant multi-mode F+10 Cepheid 
models (within the circle) in the HR diagram. The loci of other 
resonance centres are also shown. The thick and thin solid lines 
enclose the fundamental and first overtone instability strips. The 
dashed lines are the lines of the constant period (F-mode) in- 
dicated in the figure. Computations were done for set R2 of 
Table Q] 



is located between the F/IO either-or domain (hotter models) 
and a fundamental-mode-only pulsation domain (cooler mod- 
els). Therefore, from an astrophysical point of view, a double- 
periodic solution can be reached only during red-ward evolution 
by models previously pulsating in the first overtone. During the 
blue-ward evolution, models pulsate in the fundamental mode, 
which remains stable. Double-periodic pulsation is not possible. 

The properties of the computed double-periodic models, 
their periods and period ratios are collected in Table [2] (sets Rl 
and R3) and Fig. Q] (set R2). The models occupy a narrow pe- 
riod range. Compared to the observations (see Petersen diagram 
in e.g. Moskalik & Kolaczkowski 2009), the period ratios are 
slightly too low. This can be easily compensated by the decrease 
of model's metallicity. However, the main drawback of most 
of the computed models is their excessive amplitude. Only for 
models of set R3 do the amplitudes of the single-mode pulsation 
agree with the observations. What is more important, the overall 
properties of the fundamental mode and first overtone Cepheids 
are reproduced reasonably well with the convective parameters 
of set R3 (see Smolec 2009b and Baranowski et al. 2009, respec- 
tively). For the models of sets Rl and R2, the amplitudes of the 
single-mode solutions along the computed model sequences are 
always higher than observed. Particularly the amplitude of the 
fundamental mode is too high by up to 30 per cent. We will re- 
turn to the problem of excessive amplitudes later in this section. 

At first glance, the non-resonant mechanism seems to be 
the most likely cause of the double-periodic pulsation. Double- 
periodic solutions are always located between the fundamen- 
tal only and F/IO either-or domains, which is natural for the 
non-resonant mechanism. For the typical resonant excitation, the 



double-periodic domain usually emerges in the middle of the 
single-mode pulsation domain - fundamental mode domain, as 
for the 2: 1 resonance discussed by Smolec (2009a), or first over- 
tone domain, as for the three-mode resonance acting in ji Cephei 
models (Smolec & Moskalik 2007). This is because of the res- 
onant destabilisation of one of the limit cycles. Also, the non- 
resonant amplitude equations seem to capture our hydrodynami- 
cal results very well. The solid and open squares in Fig.[T]corre- 
spond to stable and unstable fixed points, computed as described 
in Sect. 12.31 through fitting the non-resonant amplitude equa- 
tions to the hydrodynamical trajectories. Except for one model 
with r e jf = 5840 K, the location and stability of the fixed points 
is consistent with the properties of the computed trajectories. 

Nevertheless, other arguments point towards a resonant ex- 
planation of the computed double-periodic models. First, all the 
computed double-periodic models are very close to the reso- 
nance centre. The proximity parameter, A (Eq. [Tol l, is given in 
Table [2] for the two relevant models of sets Rl and R3 and in 
Fig.Q~]for the models of set R2. This is not likely to happen by ac- 
cident, particularly for the three model sequences of the different 
convective parameters and the masses and luminosities. Second, 
our extensive model survey (Smolec & Moskalik 2008b) has not 
revealed any double-periodic solutions; particularly, no such so- 
lutions have been found close to the transition line between fun- 
damental mode pulsation domain and the F/IO either-or domain. 
We verified that in all models the 2d)\ = ojo + L02 resonance 
was far from the transition line. Also, the double-periodic mod- 
els were not found close to the resonance centre if the transition 
line was distant from the resonance centre. For example, models 
adopting the convective parameters of set R3 and slightly lower 
masses (4.75 M Q instead of 5.0 M ) all pulsate in the fundamen- 
tal mode when they are computed close to the resonance centre. 
This indicates that a special condition is required to generate 
the stable double-periodic pulsation - a proximity of the models 
to the resonance centre and a proximity to the discussed tran- 
sition line. Consequently, the double-periodic domain is always 
restricted to a very narrow region in the HR diagram, both in 
effective temperature and in luminosity. At the moment we can- 
not provide any explanation why these two conditions have to be 
satisfied simultaneously. 

The existence and properties of the double-periodic pulsa- 
tion depend on the values of convective parameters that enter the 
model computations. We computed several sequences of models 
with other convective parameters than R1-R3. In order to de- 
crease the model amplitudes, which are too high for sets Rl and 
R2, we computed some models with an increased eddy-viscous 
dissipation. In particular, we computed a model sequence adopt- 
ing the convective parameters of set Rl, except for the eddy- 
viscosity parameter, a m , which was set to a slightly higher value 
(a m = 0.35, instead of a m = 0.30). Unfortunately, the double- 
periodic solution simply disappeared. However, we note that it 
is hard to analyse how the existence and properties of possible 
double-periodic models depend on their convective and physi- 
cal parameters. As noted above, the double-periodic models are 
located close to the resonance centre and close to the transition 
line between the F and F/IO pulsation domains. The change in 
e.g. convective parameters while keeping the physical param- 
eters of the model fixed changes the relative location of these 
lines. Consequently, the lack of the double-periodic solutions in 
models computed close to the resonance centre may be caused 
not by the 'bad' convective parameters, but because we are in- 
vestigating a wrong part of the HR diagram. To check whether 
the inclusion of radiative losses is crucial for the existence of 
double-periodic solutions (we used y t + in all three sets with 
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double-periodic solutions), we computed two additional model 
sequences, where all parameters were left exactly the same as 
for sets Rl and R2, except for the radiative losses, which were 
turned off (y r = 0). No double-periodic models were found close 
to the resonance centre, but the transition line between the F 
and F/IO pulsation domains also shifted towards higher tem- 
peratures. Therefore, we cannot say whether the inclusion of ra- 
diative losses is necessary or not. These definite claims require a 
much more extensive model survey. Such a survey was not con- 
ducted because the double-periodic resonant solutions we found 
can model only a very small subgroup of the observed double- 
periodic F+10 Cepheids and offers no general solution to the 
reopened problem of modelling the beat Cepheid pulsation. 

Nevertheless, the presented double-periodic models are 
promising, because they are computed with a convective hy- 
drocode that includes negative buoyancy effects, and despite the 
described difficulties, their parameters agree fairly well with the 
observations. 



4. LMC double-overtone Cepheid model survey 

In contrast to our Galaxy, many more double-periodic Cepheids 
are known in the Large Magellanic Cloud. Thanks to observa- 
tional projects aimed at detecting gravitational lensing events, 
such as MA CHO (e.g. Alcock et al. [19"9"5l or OGLE (e.g. Udalski 
et al. 120081 1. a precise and extended photometry of many fields 
covering significant part of the LMC is available. The detected 
samples of Cepheids are large and homogeneous. Here we focus 
on data from the third phase of the OGLE project (Soszynski et 
al. 120081 . 

Double-overtone Cepheids are the most frequent class of 
multi-periodic Cepheids in the LMC. Two hundred and six such 
variables were found in OGLE-III data (Soszynski et al. 2008) 
compared with only 61 F+IO double-periodic Cepheids. There 
are many more single-periodic first overtone Cepheids (1238), 
but in a period range 0.5 d < Pi < 0.9 d, double-overtone 
Cepheids are the most common type of pulsators among the 
LMC Cepheids. Surprisingly, only 14 variables are confirmed 
as single-periodic second overtone pulsators. Hence, the second 
overtone is excited predominantly together with the first over- 
tone. 

The Petersen diagram for the LMC double-overtone 
Cepheids is plotted in Fig. [3] The double-overtone Cepheids 
seem to fall into two overlapping groups in this diagram. The 
first, less numerous group extends from the shortest periods up 
to periods slightly longer than P\ as 0.6 days (log P\ as -0.2). 
Within this group, one can observe a weak increasing trend in 
the P%jP\ period ratio. Members of the second, more numerous 
group that extends for longer periods (Pi > 0.6 days) form a well 
defined sequence with a decreasing period ratio. 

The evolutionary state of these objects is a matter of con- 
troversy. An analysis of the Petersen and period-luminosity di- 
agrams (Dziembowski & Smolec 2009) shows that short period 
objects are in the post-main sequence evolutionary phase, cross- 
ing the instability strip for the first time. However, objects with a 
longer period represent a challenge for the stellar evolution the- 
ory, because the inferred masses are too small for entering the 
instability strip during the core-helium burning, and the lumi- 
nosities are much higher than in the post-main sequence phase 
(see Dziembowski & Smolec 2009 , and discussion in Sect. 14. 2\ . 

A systematic survey in search for stable double-overtone 
Cepheid pulsation was not published up to date. Some 10+20 
models were computed with the Florida-Budapest hydrocode 



(see Buchler 2009), but details of these models were not pub- 
lished. Also, the exclusion of the negative buoyancy in the 
Florida-Budapest code (Sect. I2.2l i may influence the computed 
double-overtone behaviour, just as for the F+IO double-mode 
Cepheid models (Smolec & Moskalik 2008b). Therefore, it is 
very important to check whether the double-overtone models 
can be computed with the hydrocode, which includes negative 
buoyancy effects. This survey is complementary to our (unsuc- 
cessful) search for F+IO double-periodic Cepheids (Smolec & 
Moskalik 2008b), and allows for more definite statements about 
the ability to model the beat Cepheid pulsation with current con- 
vection models. The release of OGLE-III data, which contain a 
significant number of double-overtone Cepheids (Soszynski et 
al. 12008) provided additional motivation, and focused our atten- 
tion on modelling the LMC pulsators. 

4.1. Construction of models 

The modelling of double-overtone pulsation and, more gener- 
ally, pulsation involving higher order overtones, is not an easy 
task. The higher the order of the pulsation mode, the deeper 
into the envelope it penetrates. Hence, to obtain reliable peri- 
ods and period ratios, model envelopes should be deep and com- 
puted with higher resolution in the internal layers compared to 
fundamental mode Cepheid models. On the other hand, double- 
overtone Cepheids are low mass objects, characterised by small 
growth rates. This makes non-linear computations extremely 
time-consuming and therefore a relatively coarse mesh is nec- 
essary to conduct an extended non-linear model survey. In order 
to get reliable periods and period ratios with a relatively coarse 
mesh, the model structure, zoning, and depth of the envelope 
in our code were chosen to reproduce the results obtained with 
the LNA code of Dziembowski dl977l ) as closely as possible. 
Dziembowski's code, which is coupled with the Warsaw-New 
Jersey stellar evolution code (see e.g. Pamyatnykh 1999), al- 
lows the computation of deep, evolutionary models with high 
spatial resolution, providing accurate periods and period ratios 
that can be used in asteroseismic modelling (e.g. Moskalik & 
Dziembowski 2005). However, the code adopts frozen-in MLT 
convection, which is much simpler than the time-dependent 
treatment used in our code. Therefore, a comparison of the 
model periods and period ratios computed with both codes was 
conducted either for purely radiative models, or with convection 
in the frozen-in approximation in both of the codes. We note 
that with appropriate convective parameters, the KuhfuB model 
used in our envelope code can be reduced to a standard MLT 
(Wuchterl & Feuchtinger 1998 ). The resulting mesh structure for 
our envelope models is the following: the models consist of 200 
mass zones, of which 50 outer zones have an equal mass down to 
the anchor zone, in which the temperature is set to P a = 1 1000 K. 
The envelope extends down to temperatures of 8 x 10 6 K. 

In the final linear and non-linear model computations only 
one set of convective parameters was used namely set R3 of 
Table [1] The non-linear computation of double-overtone pulsa- 
tion is extremely time-consuming and prohibits a more detailed 
parameter study. Therefore, the convective parameters should be 
carefully chosen, and we believe that set R3 represents a good 
choice. We note that with the convective parameters of set R3, 
we were able to successfully model the overall properties of the 
radial velocity curves of the first overtone Cepheids (Baranowski 
et al. 120091 ). With these convective parameters, the models of the 
radial velocity curves of the fundamental mode Cepheids also 
agree with the observations (Smolec 2009b ). It is also important 
that set R3 includes the effects of radiative cooling of the convec- 
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tive elements (y t + 0). The inclusion of this effect was necessary 
to reproduce the long periods of some of the first overtone pul- 
sators (Baranowski et al. 2009}. It was also claimed that the in- 
clusion of radiative cooling was necessary to obtain the 10+20 
non-linear double-mode models with the Florida-Budapest code 
(Buchler & Kollath 2000). Finally, we would like to point out 
that the models to be discussed are hot and convection is not 
expected to be very strong. Therefore, the results should not be 
very sensitive to the exact values of the convective parameters. 

4.2. Linear modelling 

The results of the linear model survey presented in this section, 
particularly the implications for stellar evolution theory, were 
published by Dziembowski & Smolec ( [20091 DS09 in the fol- 
lowing). In the present paper, the linear models provide a back- 
ground for non-linear pulsation modelling, and are briefly sum- 
marised below. 

The pulsation models were constructed along evolutionary 
tracks computed with the Warsaw-New Jersey stellar evolution- 
ary code, which allows us to compute evolutionary phases be- 
fore core helium ignition. The computations were conducted 
for two values of metallicity that are appropriate for the LMC, 
namely Z = 0.006 and Z = 0.008. For the hydrogen abun- 
dance, X — 0.72 was adopted. The OP opacities (Seaton 2005 ) 
and Asplund et al. (2004) solar mixture were used in the opac- 
ity computations. Rotation and overshooting from the convective 
core were neglected. The mixing length parameter in the evolu- 
tionary code was set to ckmlt = 1.5. 

In addition to pulsation models computed along evolution- 
ary track (AlogL = 0.0 in the following), we have computed 
the models with artificially increased luminosity. We imposed 
higher L values (AlogL = 0.2 or AlogL = 0.4) at the bot- 
tom boundary of our models, keeping their mass fixed. The com- 
puted models cover the whole instability strip. This luminosity 
increase was intended to model either (overlarge) overshooting 
during the main sequence evolution or the core helium burning 
phase (second and third crossings of the instability strip). 

The domains of the simultaneous linear instability of the 
first and second overtone modes are plotted in the Petersen 
diagrams presented in Fig. [3] (upper panel for Z = 0.006 
and lower panel for Z = 0.008). The corresponding period- 
luminosity (PL) diagrams were published in DS09. The lin- 
ear models well reproduce the observed Petersen and PL di- 
agrams. The inferred masses of the double-overtone Cepheids 
are 3.0 ± 0.5 M . However, in order to reproduce the objects 
of longer periods (Pi > 0.6 days), which represent the major- 
ity of the sample, a significant luminosity increase compared to 
the post-main sequence evolutionary phase is required, indicat- 
ing that these objects are in the core helium burning phase. As 
noted by DS09, this represents a challenge for the stellar evolu- 
tion theory. According to current evolutionary calculations, the 
minimum mass of the star to enter the instability strip in its core 
helium burning phase is higher than 4.0 M (see references in 
DS09). For a discussion of this point we refer the reader to DS09. 

Our detailed linear analysis revealed that two types of res- 
onances can be conducive to producing stable double-overtone 
pulsation. These are 2: 1 resonance between the first overtone and 
linearly damped fifth overtone, 2cl>\ - a>s, and high-order reso- 
nance, involving the three lowest order overtones, 2a>2 = lji+ojj,. 
The loci of resonance centres within the 10+20 instability do- 
mains (and slightly beyond) are shown with asterisks (2coi = a>s) 
and pentagons (2a>2 = a>\ + 0)3) in Fig. [3] The first resonance, 
2a>i = (05, was suggested by Dziembowski (see Soszyriski et 
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Fig. 3. Domains of simultaneous linear instability of first and 
second overtones in the Petersen diagram. Black dots represent 
the observed LMC double-overtone Cepheids. Line segments are 
models with normal and artificially increased luminosity, as de- 
scribed in the bottom panel. Asterisks and pentagons mark res- 
onances between low order modes. Resonances might be help- 
ing to produce double-overtone behaviour only for the shorter- 
period Cepheids, but not for the more numerous longer-period 
ones. 



al. [2008 , DS09) as a possible factor shaping the structures visi- 
ble in Fourier decomposition parameters of the shortest period 
first overtone Cepheids. It is also well known that 2:1 reso- 
nances can be conducive to producing stable multi-mode pul- 
sation, as shown theoretically by e.g. Dziembowski & Kovacs 
(119841) . Examples of double-periodic hydrodynamical radiative 
models of both RR Lyrae and Cepheids in the proximity of a 2: 1 
resonance are also known (Kovacs & Buchler 1988 Buchler, 
Moskalik & Kovacs [19901 Smolec l2009ab . The latter resonance, 
2a>2 = 0)\ + <^3, is of the same type as the resonance operating in 
the F+IO double-periodic Cepheid models described in Sect. [3] 
We note that these resonances can be operational only in objects 
of periods shorter than Pi < 0.7 days. At longer periods, the 
non-resonant mechanism should be responsible for the observed 
double-overtone pulsation. 

4.3. Non-linear modelling 

An extensive non-linear model survey was conducted for models 
with metallicities Z = 0.006 and Z = 0.008. Model sequences 
were computed along selected sequences described in the previ- 
ous section chosen to cover the significant part of the Petersen 
diagram. 

For metallicity Z = 0.006, six non-linear model sequences 
were computed. Three sequences are located within the instabil- 
ity domains along evolutionary tracks (A log L = 0.0) for stars 
of masses 2.5 M , 3.0 M and 3.5 M (solid lines in Fig. |3}. 
The other three are located along horizontal paths with arti- 
ficially increased luminosity (2.5 M Q , AlogL = 0.4; 3.0 M , 
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Fig. 4. Petersen diagram showing the location and stability in- 
formation for the computed non-linear convective models. The 
models were chosen to cover the significant part of the Petersen 
diagram. These are models of 2.5 M , 3.0 M and 3.5 M located 
along evolutionary tracks (A log L - 0.0) and models with arti- 
ficially increased luminosity of masses 3.5 M Q (AlogL = 0.2), 
and 2.5 M and 3.0 M (AlogL = 0.4). All computations are 
for the metallicity parameter Z = 0.006. Black dots represent 
the observed LMC double-overtone Cepheids. The only double- 
periodic models are found near the 2u)\ = a>s resonance centre, 
which is marked with an arrow. 
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Fig. 5. Same as Fig. |4] but for metallicity parameter Z = 0.008. 
The chosen models cover the significant part of the instability 
strip. These are models of 2.5 M , 3.0 M and 3.5 M located 
along evolutionary tracks (A log L - 0.0) and models with an 
artificially increased luminosity of masses 2.5 M and 3.0 M 
(A log L = 0.4). The only double-periodic models are found near 
the 2ll>i = u>5 resonance centre, which is marked with an arrow. 



AlogL = 0.4; 3.5 M , AlogL = 0.2). The location of all the 
computed models is shown in Fig. [4] 

Five non-linear model sequences were computed for a higher 
metallicity, Z = 0.008. Three sequences correspond to first 
crossing evolutionary models of masses, 2.5 M , 3.0 M and 
3.5 M . The remaining two sequences correspond to horizontal 
paths with an artificially increased luminosity (3.0 M , A log L = 
0.4 and 3.5 M , AlogL = 0.4). The location of the computed 
models is shown in Fig. [5] 

For each model in a sequence five hydrodynamical integra- 
tions were conducted, each initialized with a different mixture of 



the first and second overtone's velocity eigenvectors. For most of 
the models the integrations were carried over 8000 pulsation cy- 
cles. For the least massive models located along the evolutionary 
track (2.5 Mg, AlogL = 0.0), integrations were twice as long, 
because the growth rates are very low for these models. Some 
individual trajectories for models located along the evolutionary 
track with 2.5 M (Z = 0.006) are presented in Fig. [6] 

The modal selection information for the computed models is 
presented with different symbols in Figs. [4] and [5] It was derived 
through the analysis of the computed trajectories (presented e.g. 
in Fig. [6j. As described in Sect. I2.3I hydrodynamical computa- 
tions were not analysed with amplitude equations because of the 
resonances. We also note that for the discussed models the fun- 
damental and third overtone modes can be linearly unstable as 
well. We checked with the analytical signal method that these 
two modes are present during the transient evolution, but they 
decay, and except for one model (see next paragraph) are not 
present in full amplitude pulsation. 

The results are qualitatively the same for both metallicities, 
and below we focus our discussion on the lower metallicity mod- 
els (Figs. |4] and [6j. The most interesting modal selection sce- 
nario is visible along the 2.5 M evolutionary track. The individ- 
ual trajectories of some models in this sequence are displayed 
in Fig. [6] For the hottest model in this sequence (not shown 
in Fig. [6j, 20/30 hysteresis is possible. Depending on the ini- 
tial conditions, the trajectories evolve either towards the second 
overtone attractor or towards the third overtone attractor. In six 
consecutive cooler models, only single-periodic pulsation in the 
second overtone is possible. Then, a double-overtone domain 
emerges. A multi-mode attractor is clearly visible for models 
with logr eff = 3.8349, 3.8316, and 3.8284. Periods and pe- 
riod ratios for these models are given in Fig. [6] The only model 
redward of this domain pulsates in the first overtone. Thus, the 
double-overtone domain appears in between the first overtone 
and the second overtone pulsation domains. Although this is a 
typical scenario for a non-resonant mechanism, a resonant mech- 
anism cannot be excluded a priori. The arrow in Fig. |4] shows 
the location of the 2a>i - a>s resonance centre. It is located 
quite close to the double-overtone domain. The resonant desta- 
bilisation of the first overtone limit cycle, if it occurs at a proper 
position along the model sequence, can in principle lead to the 
emergence of the double-overtone domain instead of the either- 
or 10/20 domain. Note that the 2:1 resonance can affect the 
pulsations in a very wide range of period ratios, as is observed 
for bump Cepheids, and that destabilisation may not necessarily 
occur exactly at the resonance centre (see Smolec 2009a). We 
postpone the discussion of the nature of the computed double- 
periodic models to the next paragraphs in this section. 

The remaining five non-linear model sequences of Z = 0.006 
all display qualitatively the same mode selection scenario, which 
is different however from scenario described for the 2.5 M se- 
quence. A double-periodic domain is not present. Instead, the 
either-or 10/20 domain is located in between the single-periodic 
first overtone and second overtone pulsation domains. The dis- 
cussed resonances seem to have no visible effect on the modal 
selection along these sequences. 

For higher metallicity models, the overall modal selection is 
very similar to the lower metallicity models (compare Figs.|4]and 
|5j. Again, the most interesting results are obtained for models of 
the shortest periods with masses of 2.5 M . Here, a multi-mode 
domain emerges. Simultaneous pulsation in the first and second 
overtones is possible for three models along the sequence. The 
centre of the 2:1 resonance between first and fifth overtones, 
marked with arrow in Fig. [5] falls exactly in the middle of the 
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Fig. 6. Fractional radius amplitude for the second overtone (A2) and for the first overtone mode (A,), for nine consecutive models 
located along evolutionary track of 2.5 M Q (Z = 0.006). The models were initialized with a range of amplitude ratios. The effective 
temperatures of the models and linear P5/P1 period ratios are given in each panel. The approximate location of the attractors is 
marked with filled circles. For the models with a double-periodic attractor the periods P, and period ratios, P2IP1, are also given. 
All trajectories in the first five panels evolve towards a single-periodic second overtone attractor. In the next three panels (panels 
6-8) a double-periodic attractor is clearly visible. In the last panel, all trajectories evolve towards a single-periodic first overtone 
attractor. 



double-periodic domain. This supports the hypothesis that the 
resonance is crucial in establishing the stable double-periodic 
behaviour that we found for both metallicities. For the remain- 
ing four model sequences with metallicity Z = 0.008, no traces 
of multi-mode behaviour are found. Again, as for the Z = 0.006 
models, it seems that resonances do not affect the modal selec- 
tion along model sequences of higher masses (longer periods). 

Considering the double-overtone models we found it is hard 
to judge which mechanism, resonant or non-resonant, is respon- 
sible for the computed beat pulsation. Both double-periodic do- 
mains (for Z = 0.006 and Z = 0.008) are located very close 
to the 2<x>\ = cl>s resonance centre. For the Z = 0.008 models 
the resonance centre falls exactly in the middle of the double- 
periodic domain. This suggests the resonant origin of the com- 
puted beat pulsation. On the other hand, these double-periodic 
domains are located in between two single-periodic pulsation 
domains, first overtone pulsation domain (to the red) and sec- 
ond overtone pulsation domain (to the blue). This is typical for 
the non-resonant scenario. Also, for all studied sequences of 
masses M > 2.5 M , which cross the 2w, = 0)5 resonance cen- 



tre (see Fig.O, no traces of beat pulsation were found. Instead, 
an 10/20 either-or domain is present. Although the presence of 
the resonance may be not sufficient to excite the beat pulsation 
and other factors may be necessary (like for the F+IO models 
discussed in the previous section), it is hard to identify these fac- 
tors and prove their necessity. Therefore, we do not presume to 
know which mechanism, resonant or non-resonant, underlies the 
double-overtone behaviour we found. 

Although some double-overtone models were found, the 
overall results are not satisfactory. The computed double- 
overtone domains are narrow and are located at P, w 0.3 d, i.e. 
at the short-period end of the observed domain of the double- 
overtone pulsation in the LMC. No double-overtone models 
were found at longer periods, particularly in a period range, 
0.5 d < Pi < 0.9 d, where the double-overtone pulsation is the 
most common form of pulsation in the LMC. Surprisingly, in 
many models only pulsation in the second overtone is possible. 
Observationally, this form of pulsation is very rare, as only 14 
such objects are identified in the LMC (Soszyhski et al. 2008 
still, observational selection can be a factor here). Therefore, 
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also for the 10+20 double-overtone Cepheids the convective 
hydrocode fails to reproduce the observed modal selection. This 
conclusion is not as strict as for the F+IO Cepheids, because 
only one set of convective parameters was explored. On the other 
hand, for hot overtone models, convection is not expected to play 
a crucial role, and results should not depend strongly on the val- 
ues of convective parameters and/or effects included in the con- 
vective model. 

5. Conclusions 

The longstanding problem of non-linear modelling of beat 
Cepheid pulsation remains open. Our search for stable double- 
periodic Cepheid pulsation with convective hydrocode, which 
correctly includes negative buoyancy effects, yielded unsatisfac- 
tory results. This search is a difficult and time-consuming task. 
The turbulent convection model we use contains several free pa- 
rameters. Also, the Cepheid domain extends over a wide range 
of luminosities in the HR diagram, covering different masses and 
evolutionary stages. Consequently, extensive model computa- 
tions are needed. Our recent model surveys (Smolec & Moskalik 
2008b, Smolec 2009b, this paper) cover different stellar systems 
(Galaxy, LMC) and different evolutionary stages (post-main se- 
quence and helium burning objects). Different sets of convective 
parameters were investigated. An extensive parameter study was 
done for F+IO Cepheid models (Smolec & Moskalik 2008b). 
For the 10+20 models presented here, only a limited study with 
carefully chosen parameters was possible, owing to the longer 
computation time needed for these models (see Sect. 4.1). Only 
a few double-periodic models, limited to a narrow period ranges, 
were found. In particular, we have not found any large domain 
in which the non-resonant mechanism is the sole cause of the 
double-mode pulsation. 

In most of the observed Cepheids a non-resonant mecha- 
nism should be responsible for the simultaneous pulsation in two 
modes. With purely radiative hydrocodes, non-resonant double- 
mode models are found only for RR Lyrae stars, and only if the 
artificial viscosity is significantly reduced (Kovacs & Buchler 
1993). These models are however sensitive to numerical de- 
tails and their pulsation amplitudes are too high. On the other 
hand, several computed resonant radiative beat Cepheid and 
RR Lyrae models (Kovacs & Buchler[1988 Buchler, Moskalik 
& Kovacs 1990, Smolec 2009a) disagree with the observations. 
The inclusion of turbulent convection into pulsation hydrocodes 
seemed to offer the solution. Using modified KuhfuB convec- 
tion model where negative buoyancy effects were excluded (e.g. 
Kollath et al. 1998 ), the Florida-Budapest group published many 
non-resonant beat Cepheid models (e.g. Kollath et al. 2002, 
Buchler 2009). However, the exclusion of the negative buoyancy 
is physically not justified and, as shown by Smolec & Moskalik 
(2008b]), it is the main cause of the computed double-mode pul- 
sation. Using the original KuhfuB prescription including the neg- 
ative buoyancy, we were unable to find any non-resonant F+IO 
double-mode Cepheid models (Smolec & Moskalik 2008^). In 
this paper we conducted an additional search for double-periodic 
Cepheids, including double-overtone (10+20) pulsators. With 
our hydrocode, which includes negative buoyancy effects, we 
found several double-periodic models of both F+ lO and 10+20 
type. In the first case, the double-periodic pulsation is most 
likely caused by the 2<x>\ = ojq + toi resonance. In the latter 
case, the exact mechanism cannot be identified beyond doubt. 
We only note that the 2:1 resonance, 2co\ = co^ may be opera- 
tional in these models. The computed models are very interest- 
ing. Their periods and period ratios agree well with the observa- 



tions. Nevertheless, they are restricted to very narrow domains 
in the HR diagram and offer no general solution to the problem 
of modelling the beat Cepheid pulsation. 

Turbulent convection is an important phenomenon deter- 
mining the properties of Cepheid envelopes and cannot be ne- 
glected in pulsation models. However, the numerical expense 
of non-linear computations requires relatively simple convective 
recipes. So far, only one-equation models for the generation of 
turbulent convection were implemented in pulsation hydrocodes. 
Two such models were used in the computation of Cepheid 
models, the KuhfuB (1986) model and the Stellingwerf model 
(Stellingwerf [19821 Bono & Stellin gwerf [19921 . O ur exten sive 
computations (Smolec & Moskalik l2008bl Smolec |2009bl this 
paper) indicate that the original KuhfuB model (i.e. including 
negative buoyancy) is incapable of reproducing the majority of 
the observed beat Cepheids. In the original Stellingwerf ( 1982) 
model, the negative buoyancy effects were neglected, because 
in that model the turbulent source function is proportional to 
square-root of the superadiabatic gradient, S oc y/Y. This makes 
the original Stellingwerf model very similar to the modified 
KuhfuB model used in the Florida-Budapest hydrocode, and in- 
deed, as investigated by Buchler & Kollath (2000) both convec- 
tion recipes lead to qualitatively the same results. Stellingwerf's 
(119821 ) functional form of the source function was criticised by 
Gehmeyr & Winkler (U992| > who pointed out the arising prob- 
lems associated with the long time-scale of decay of the tur- 
bulent eddies in regions that became convectively stable dur- 
ing pulsation. Bono & Stellingwerf (1992) modified the origi- 
nal Stellingwerf recipe through setting S oc sgn(F) y/\Y~\- With 
this model for negative buoyancy, plenty of Cepheid mod- 
els were computed with the Italian code (e.g. Bono, Marconi 
& Stellingwerf I2000T ). However, no double-periodic classical 
Cepheid models were ever found (Bono, private communica- 
tion). This supports our results. Including the negative buoyancy 
in the convective models (both in KuhfuB and in Stellingwerf 
models), wide domains of double-periodic pulsation cannot be 
found. We conclude that with the one-equation models used in 
non-linear Cepheid modelling, double-periodic Cepheid pulsa- 
tion cannot be modelled satisfactorily. 

The question of how the hydrocodes should be modified to 
solve the puzzle of beat Cepheid pulsation is difficult to an- 
swer. A simple convection model seems the most severe short- 
coming of the present hydrocodes. More sophisticated convec- 
tive recipes are being developed (Stokl 2008 Buchler 2009), 
however, many problems have to be solved before they can 
be applied in non-linear modelling. Another shortcoming of 
the present models is the simple structure of the constructed 
Cepheid envelopes - they are chemically homogeneous and non- 
rotating. Without detailed model computations however, it is 
hard to judge whether e.g. possible differential rotation could 
have an effect on modal selection. The treatment of radiation is 
also very simple in most of the hydrocodes, i.e. diffusion ap- 
proximation is used (our code, the Florida-Budapest code and 
Italian code). However, a more detailed time-dependent treat- 
ment did not lead to substantial changes in the computed models 
(Feuchtinger, Buchler & Kollath 2000). Finally, in the present 
models only radial pulsation is considered. Non-radial modes are 
neglected. There is a growing evidence that non-radial modes are 
excited in classical Cepheids, including the beat Cepheids (see 
e.g. Moskalik & Kolaczkowski 2009). The coupling between 
non-radial and radial modes can affect the modal selection, and 
in our opinion, this is one of the ideas that should be investigated 
first. For this purpose, an analytical approach based on the am- 
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plitude equation formalism may be used, without the necessity 
of developing a non-linear non-radial hydrocode. 
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